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Abstract 


We address the question which additional information on the source shape 
and dynamics can be extracted from three-particle Bose-Einstein correlations. 

For chaotic sources the true three-particle correlation term is shown to be sen¬ 
sitive to the momentum dependence of the saddle point of the source and to its 
asymmetries around that point. For partially coherent sources the three-pion 
correlator allows to measure the degree of coherence without contamination 
from resonance decays. We derive the most general Gaussian parametriza- 
tion of the two- and three-particle correlator for this case and discuss the 
space-time interpretation of the corresponding parameters. 

PACS numbers: 25.75-|-r, 13.85 Hd, 24.10-i. 

I. INTRODUCTION 

Two-particle Bose-Einstein interferometry (also known as Hanbury Brown-Twiss inten¬ 
sity interferometry) as a method for obtaining information on the space-time geometry and 
dynamics of relativistic heavy ion collisions has recently received intensive theoretical and 
experimental attention. Detailed theoretical investigations (for a recent review see Ref. [Q) 
have shown that high-quality two-particle correlation data can reveal not only the geometric 
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extension of the particle-emitting source but also its dynamical state at particle freeze-out. 
This information is encoded in the second central space-time moments of the “emission 
function” S{x, K), i.e. of the Wigner phase-space density of the source. For chaotic sources, 
certain linear combinations of these moments can be extracted from the two-particle cor¬ 
relation function C 2 {q,K) by htting it to a Gaussian in the relative momentum q of the 
pair H- These second space-time moments give the size of the regions of homogeneity 
which effectively contribute to the emission of particle pairs with a given pair momen¬ 
tum K] collective dynamics of the source results in a characteristic iF-dependence of these 
homogeneity regions 

More detailed information on the space-time structure of the source may be hidden 
in possible non-Gaussian features of the correlation function C 2 {q, K) even if they are hard 
to extract; due to the symmetry under q —> —g, however, only even space-time moments of 
the source are accessible via two-particle correlations. In this paper we will extend previous 
studies of multi-particle correlations P-p!^ and show that three-pion correlations provide 
in principle additional information on the space-time characteristics of the source which 
cannot be obtained from two-particle interferometry. We show in particular in Sec. ^ that 
for completely chaotic sources the true three-pion correlations are determined by the phase 
of the two-particle exchange amplitude [^,0 which drops out from the two-particle cross 
section. This phase is shown to be sensitive to the rate at which the saddle point x{K) of 
the source, from which most pairs with momentum K are emitted, moves as K changes, and 
to the asymmetries of the emission function around this saddle point via its third central 
space-time moments. Unfortunately, this phase turns out to be generically small, and its 
sensitivity to these asymmetries is very weak, making them extremely hard to measure. 

In the absence of such a non-trivial phase, three-particle correlations can still be used 
to test the chaoticity of the emitting source. To this end we derive the expressions for 
two- and three-particle correlations for chaotic and partially coherent sources and establish 
their respective relationships. Our treatment differs from previous studies of multi-particle 
Bose-Einstein correlations in that we consistently express the correlation functions through 
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the source Wigner density, even for partially coherent sources. This enables us to relate the 
shape of the correlators as functions of the various relative momenta to certain space-time 
features of the source. To the best of our knowledge the corresponding relations for partially 
coherent sources (Eqs. (^) - (|36D) are new. 


II. CHAOTIC SOURCES 


For a chaotic source, the two-pion correlation function C 2 {pi,Pj) can be expressed as 

Bin 

P2{Pi,Pj) 


C2{Pi,Pj) = 


Pl{Pi) Pl{Pj) 


= 1 + 


1/ d^x S(x, Kij) 


= 1 + 


\Pi. 


( 1 ) 


fd^xS(x,Pi) fd^yS(y,pj) pupjj 

Here P 2 {pi,Pj) is the two-pion inclusive cross section, and Pi{pi) is the single-particle inclu¬ 
sive spectrum. S{x,p) is the single-particle Wigner density of the source, i.e. the quantum 
mechanical analogue of its phase-space distribution. The average and relative 4-momenta 
Kij = {Pi + Pj)/‘^ and qtj = Pi — Pj satisfy the constraint qij ■ Kij = 0 which results from 
the on-shell nature of the observed momenta pi. The two-particle exchange amplitude pij is 
dehned as [ITBIP^ 


pij = PiQij^Kij) = ^EiEj{d\pi)d{pj)) 

= [ d^xS{x,Ki,)E^--^ = fi,E^-. 


( 2 ) 


From (0) it follows that pij = p*i and thus fij = fji and (j)ij = —(j)ji. Correspondingly, 
(t>n = 0, Pa = fa, and fij must be an even function of qij while fij is odd in qij. 

The single-pion spectrum can be written as 


Pi{Pi) = J d^^xSix^pi) = fi, 


( 3 ) 


while the true two-pion correlation function is dehned by 


R2{i,j) = R2{Pi,Pj) = C2{pi,Pj) - 1 = 


J 2 




fiif' 


( 4 ) 


jj 
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Similarly, the true three-pion correlation function is given by ||l3|,p!6|-pl 


R3(,Pi,P 2,P3) = C3{pi,P2,P3) - ^ 2 ( 1 , 2 ) - i?2(2,3) - i?2(3, 1) - 1 
Re (pi 2 P 23 p3i) 


= 2 
= 2 


fii /22 /as 
/12 /23 /3I 


( 5 ) 


COS(012 + 023 + 03l) • 


/ll 022 033 

Since the real parts fij of the exchange amplitudes pij can be extracted from the two- 
pion correlator, for chaotic sources the only additional information contained in the 3-pion 
correlation function resides in the phase |]T^ 


•h = 012 + 023 + 031 ; 


( 6 ) 


it is a linear combination of the phases of the three exchange amplitudes pi 2 , P 23 , and P 31 
which enter the true 3-pion correlator R 3 . This phase is odd under interchange of any two 
particles. It can be isolated by normalizing R 3 with respect to the true 2-pion correlator R 2 : 

R3{Pi,P2,P3) 


r3{PuP2,P3) = 


= 2 cos $ . 


(7) 


^R2(1,2)R2(2,3)R2(3,1) 

In order to understand which space-time features of the source affect the phase $ (and 
thus the normalized true 3-pion correlation function r^) we expand the exchange amplitude 
Pij for small values of qij = Pi — Pj . We dehne the average of an arbitrary space-time 
function f{x) with the source distribution S{x,Kij) as 

Jd'^xfix) S{x,Kij'^ 


Jd^xS{x,Kij) ■ 

This average is a function of the pair momentum Kij. Using (|^) we thus get 
Pij = Pi{Kij) 1 + i{qij-x)ij - ^ (^{qij-xf),, - ^ (^{qij-xf),, + O (gj) 
Separating real and imaginary parts we hnd, after a little algebra. 


fij = Pi{Kij) 1 - ^ (^{qij-Xij)^),, + O (gj) 


( 8 ) 


(9) 


( 10 ) 


and 
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where 


4)ij = qij-{x)ij - ^ + O (q^^ , 


( 11 ) 


Xij = X - {x)ij = X - x{Kij) 


( 12 ) 


is the distance to the “saddle point” of the sonrce, i.e. to the point of maximnm emission 
for pions with momentnm Kij. According to Eqs. ([T^) and (§), the two-pion correlator 
is sensitive to the second central (i.e. saddle-point corrected) space-time moments of the 
emission fnnction S{x,Kij) [@,||, with higher order corrections from all even central space- 
time moments. The phase <h, on the other hand, contains information on the odd space-time 
moments. Expanding S{x,Kij) aronnd the average momentnm K of the pion triplet, 

P1+P2+ P 3 Ki2 + K23 + A'gi 


K = 


1 


K-ij — K ^ (jlik + qjk) 1 * 7^ J 7^ ^ ) 

and nsing qi 2 + q '23 + qzi = 0, we hnd from Eqs. (||) and (|lTl) 


(13) 

(14) 


2 ^ 1 " ?23 


d{x^) d{x^) 


dK^ 


dK^ 


24 L 


U U A \ U IP A 

(Il2<ll2<l23 + 923923^12 


(x^) d^jxu) d^{xx) 


dK’^dK^ dK^dKi^ dRf^dK’^ 


(15) 


2?12?23(912 + fe)^ {x^XyXx) + 0{q^) . 


Here the average withont snbscripts 


U{x)) = 


Jd^xf{x)S{x,K) 


(16) 


/ d^x S{x, K) 

denotes the space-time average with the emission fnnction evalnated at the mean momentnm 
K of the pion triplet, and 


X = X — {x) = X — x{K). (17) 

Eq. (^) is the main new resnlt of this Section. One easily checks that has it the correct 
symmetries nnder particle exchange. It shonld be noted that, dne to the on-shell constraint 
qij ■ Kij = 0, only three of the fonr components gh are independent. The resnlting relation 
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(? Qij ' f^iji with. l3ij KijI(^K ^ij , 


( 18 ) 


can be used to eliminate the redundant g-components in Eq. (p!5D, thereby mixing spatial 
and temporal components of the corresponding coefficients. This is a well-known problem 
also for the two-pion correlator (see, e.g., |^) which prohibits a clean model-independent 
separation of the spatial and temporal widths of the source. 


Eq. (pffiD features two types of contributions to the phase <h: The formally leading con¬ 
tribution enters at second order in the relative momenta qij and is proportional to the rate 
dXfj,{K)/OK'^ with which the saddle point of the emission function changes as a function of 
the pion momentum K. This term will in general be non-zero even for emission functions 
with a purely Gaussian x-dependence. It gives rise to a leading g^-dependence of the normal¬ 
ized true three-particle correlator = 2 cos <h. At order g^ the phase <I> receives additional 
contributions from the second A'-derivatives of the saddle point as well as from the third 
central space-time moments {x^XyX\) of the source. The latter are the leading contributions 
from a possible asymmetry of the emission function S{x,K) around its saddle point x{K)] 
they vanish for purely Gaussian emission functions. We see that they enter the normalized 
three-particle correlator at order g^ in a mixture with the A'-dependence of the saddle 
point. This renders their isolation essentially impossible. 

In contrast to the widths of the emission function, which affect the two-pion correlator 
at second order in the relative momentum, the additional structural information which can 
(in principle) be extracted from the (normalized) three-pion correlator is seen to enter at 
most at fourth order in g. Their measurement is thus very sensitive to an accurate removal 
of all leading g^-dependences by proper normalization to the two-particle correlators. To 
achieve this looks like an extremely difficult experimental task. We are therefore somewhat 
pessimistic about the short-term prospects of extracting additional structural information 
about the source from three-pion correlations. 

If the phase $ and the information it contains about the source are inaccessible, what 
else can three-pion correlations be used for experimentally? The answer is that one can test 
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the assumption that the source is chaotic. This has been pointed out previously in Refs. 

where specihc simple parametrizations for the two- and three-particle correlators (as 
well as for higher order correlations) were assumed and the relationship between the various 
parameters was studied. We will here derive more general expressions which, in principle, 
permit such a test without making any simplifying assumptions about the shape of the 
source. 

Before proceeding to the discussion of Bose-Einstein correlations from partially coherent 
sources, we would like to close this Section with a few short remarks on the effects from 
resonance decays. It is well known that partial coherence in the source leads to 

incomplete correlations in the two-particle sector, in the sense that R 2 {q,K) at vanishing 
relative momentum q = 0 does not approach the ideal value R 2 { 0 ,K) = 2 for chaotic 
sources. In actual experiments there are, however, other possible reasons for apparently 
incomplete two-particle correlations. Most importantly, pions from the decay of long-lived 
resonances contribute to the correlator only at very small values of q and thus (due to limited 
2 -track resolution) may escape detection in the correlation signal while fully contributing to 
the single-particle spectrum, thereby reducing the apparent correlation strength even for a 
completely chaotic source [P^-|27[|. In a Gaussian parametrization of the exchange amplitude 
this can be implemented by writing instead of Eq. (|T^) for qij ^ 0 


11,19 


/p- 


P,{K,,) exp 





(19) 


where, up to second order in q, Rf^'^{Kij) = with the source average on the r.h.s. 

being taken only over the “core” of pions from direct emission and from the decays of 
short-lived resonances The two-particle correlator then becomes 


^2(h j) = HKij) 


exp 




( 20 ) 


PM P,ip,) 

and for vanishing relative momenta q the three-particle correlation function assumes the 


value 


C 3 {pi=P 2 =P 3 =K) = 1 + 3 X{K) + 2 XP^{K ). 


( 21 ) 
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Note, however, that the expression (|^ for the normalized true three-pion correlation function 
is not affected by resonance decay contributions and remains unchanged. This will no longer 
be true for partially coherent sources. 


III. PARTIALLY COHERENT SOURCES 


Expressions for the n-particle inclusive spectra from partially coherent sources have been 
previously derived, with differing methods, in Refs. [|I^|l6H^. In the covariant current 
formalism of Refs. one decomposes the classical source current which creates the free 

pions in the hnal state into a coherent and a chaotic term: 


T(x) Jcohi^^') T '-^cha(^) • 


( 22 ) 


Following the treatment of Ref. this leads to the following dehnition of the single-particle 
Wigner density (“emission function”) of the source: 


S(i, K) = Ie-'^-y(r(x + l)J(x-%')) 


(23) 


with 


So„h(x, K)=j e-“y-> J:Jx + f) j„h(i - I), 


Sa,(x,K) = 


2(27r) 

^ y a-iK-y / j* 


2(27r) 


(^cU^ + i)^cha(a;-f)). 


(24a) 

(24b) 


The average on the r.h.s. of the dehnition (P4b|) for the chaotic part of the emission function 


is dehned as in Ref. [15|, and we used 


(^cha(a^) -Jcohiy)) = 0. 


(25) 


The Wigner density of the full source is thus the sum of a coherent and a chaotic contribution; 
no mixed terms occur because the chaotic and coherent source currents do not interfere. 
This allows to carry over the intuitive and very successful Wigner function language for 
fully chaotic sources to the case of partially or completely coherent sources. 
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We now write 


A, = d^xS{x,Ki^)e^'i^^-^ 


— I coh _ p . pi^ij I f, , pi<t>ij 

rij ' Fij — -f ij c W c , 


(26) 


where Kij ^ (p; + Pj)/2, (jij =Pi- pp and 




(27a) 

(27b) 


As shown in Ref. this yields the two-pion correlation function in the form 

+ 2fijFij cos{^ij - (pij) 


C2iPi,Pj) = 1 + R2ihj) = 1 + 


while the three-particle correlation is given by 


if a + Rii)ifjj + Fjj) 


(28) 


C3iPi,P2,P3) — 777— TFF7~~F ^/—7 ~ 1 + -^ 2 ( 1 , 2) -|- i?2(2, 3) -|- i?2(3,1) 


PiiPi) PliP2) PiiP3) 
2 


PliPl) PliP2) PliP3) 

+ /l2-^23-^31 COS(012 -|- $23 + $3l) 

+ Fi2f23p31 COs($i2 -|- 023 + $3l) 

+ Fi2F23f3l COs($i2 -|- $23 + 03l)) • 


{F12F23F31 COs($i 2 + $23 + $31) 


(29) 


Similar expressions were derived in Ref. |T^. The two- and three-particle correlations are 
seen to vanish for completely coherent sources (Fjj —>• 0 Vi,j). In the opposite limit (fij —>• 
0 Vi, j) one recovers the results from Sec. || for completely chaotic sources. 

The representations (|^) and (|^) permit us to write down for and $jj, 0*^ similar 

small-g expansions as in Eqs. (0) and (pTl) ; the corresponding averages are dehned with 
respect to the chaotic and coherent parts, respectively, of the Wigner function (^). In 
the true two-pion correlation function R 2 {i,j) of Eq. (^^, the hrst term thus contains 
information on the second central space-time moments of S'cha((^, Kij) while the second term 
mixes the second moments of Scha.ix, Kij) and Scohix, Kij) in a rather nontrival way. Since 
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the number of measurable parameters in R 2 {i,j) is the same as before, this implies a relative 
loss of information: the second space-time moments of Scha and S'coh can neither be separated 
nor do they simply combine to the second central moments of the total source S = Sr ha +Smb- 
This complication goes hand in hand with a similar one in the three-pion correlator: 
Dehning the true three-pion correlator as before, 


i?3(l, 2, 3) = Cs{pi,p2,P3) - 1 - 2) - i?2(2, 3) - i?2(3, 1) 

2 


(/ll + Fii){f22 + F 22 ){f 33 + F33) 

X {F 12 F 23 F 31 COs(<hi2 + $23 + $ 31 ) + /l2-^23-p31 COs(0i2 -|- $23 + $3l) 

+ -^" 12 / 23-^31 C0S($12 + 023 + $ 31 ) + -^ 12 -^ 23/31 C0S($12 -|- $23 + 03l)) , (30) 

one sees that, in contrast to Eq. (^) for chaotic sources, the phase factors can no longer 
be isolated by normalizing R3 with a proper combination of two-particle correlators i? 2 - 
This means that, in a samll-g expansion, i?3(l,2,3) contains leading terms of second order 
in q which are independent of those occurring in the two-particle correlator. On the one 
hand, those terms supplement the incomplete information from R2 on the second space- 
time moments of the source; on the other hand, they render the measurements of source 
asymmetries impossible. 

The full reconstruction of all the (in principle) measurable information obviously requires 
a measurement of i? 2 (b j) and i? 3 (l, 2, 3) as a function of all nine components of Pi,P 2 ,P 3 ■ 
In view of the technical complexity (both experimental and theoretical) of such a program 
this is not likely to happen soon. It must, however, be mentioned that simple one- or two- 
parameter Gaussian parametrizations as suggested in Refs. ||I7|JI^, [I^ are not sufficient for 
this purpose because they very strongly prejudice the form of the source. 

To pursue this last point a little further, let us define the (momentum-dependent) chaotic 
fraction of the single particle spectrum 

Fa J Scha{x, Pi) 


^iPi) = 


fii + Fii Jd'^xS{x,pi) 


(31) 
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The coherent fraction is accordingly fa/{fa + Fa) = 1 — e{pi). For vanishing relative mo¬ 
mentum qij = 0 {i,j = 1, 2, 3), we then have 


R 2 {p,p) = e{p){2 - e(p)), 

^ 3 (P,P,P) = 2e^(p) (3 - 2e(p)). (32) 


For completely chaotic sources, e(p) = 1, we recover the results of Sec. |T|. For partially 
coherent sources, the normalized three pion correlator at vanishing q is given by 


, ^ R3{p,P,p) o FTT (3 - 2e(p)) 


(33) 


which, in general, deviates from the chaotic limit T^{p,p,p) = 2. 

It would thus seem to be a simple matter to check the limits of R 2 and i ?3 for vanishing 
relative momenta and construct the ratio (|S3| ) in order to see whether or not the source 
contains a coherent component. In practice, however, the g = 0 limit can not be measured 
directly, but requires an extrapolation of data at hnite q to zero relative momenta. It is well 
known that such an extroplation can be very sensitive to the assumed functional behavior 
of the correlator at small q. As we will now show our results provide a basis for a reasonable 
parametrization of R 2 and R^ for small q. 

To this end we start from Eqs. (pSD and (|30|) together with the small q expansions (pTl), 
(0). Noting that R 2 must vanish for g —> 00 , a parametrization which is correct up to 
second order in g is given by 

R 2 {i,j) ~ exp [-q^jq^jRf,^{Kij) 

+ 2e{Kij){l - e{Kij)) exp [-|g-)g)'j {R/^uiFij) + r^^{Kij))\ cos {qij-s{Kij)) . (34) 


It follows from Eqs. (|l0|) , (|TI[) that here 



_ /yM z. 1 ^ \cha 
/ ij ? 

(35a) 


/ IJ ? 

(35b) 


= . 

(35c) 
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Eq. (^) neglects an additional factor P‘^{Kij)/P{pi)P{pj) which is unity for exponentical 
single particle spectra 0. Eq. (|3^ ) differs from the parametrization suggested in Ref. 


by the factor cos {qij ■ s{K)) exp 




~tJ.v{Kij) ; the parametrization of Ref. [|^ is thus 
not general enough. (It essentially assumes that the coherent part of the source is pointlike 
(in space and time!) and localized at the saddle point of the chaotic part of the source.) 
Note that from Eq. (^) one must still eliminate the redundant g-component via the on-shell 
constraint (|T^). 

The three-pion correlator can similarly parametrized as 

R 3 {Pi,P 2 ,P 3 ) = 2 e^{K) exp [- (^^ 2^2 + ^ 23^23 + ^^ 12^23 + Rt^v{K) 


X 


e{K) 

+ (1 - e{K)) cos(gi2-s(R:)) exp [\ql2ql2{R^lu{K) - 
+ (1 - e{K)) cos{q23-s{K)) exp (|g 23 ? 23 (^Mi^(^) “ 

-f (1 — e{K)) cos((gi2 + g3i)-s(R')) 

X exp (i((/i 2 + <l2s)'‘(<li2 + ‘l23)‘'{RpAK) “ . 


(36) 


This again generalizes the parametrizations given in Refs. |]^^|; according to Eqs. (PUD, 
(|n|), it is correct up to the second order in q if one approximates Pf{Kij/Pi{pi)/Pi{pj) ^ 1 


as well as e{Kij) ^ ^{K)- The parametrizations of Ref. are recovered in the limit 

of a pointlike coherent source, r^y{K) = 0, and assuming s{K) = 0. (The first of these 
two assumptions is explicity stated in Ref [0.) One can easily convince oneself that at 
qi2 = 0, for example, the term cos{q23-s{K)) exp \q 23 q 23 {Rij.v{K) - 'r^ix(R')) enters R 3 {q 23 ) 
with a different weight than R 2 {q 23 )- Thus R 3 provides additional information which allows 
to separate R^y{K) from r^y{K) and thereby the widths of the chaotic and coherent parts 
of the source. 

In practice, one must also take into account resonance decays. Since it follows from the 
discussion at the end of in Sec. that the longlived resonances do not affect the intercept 
(p3|) of the normalized true three-pion correlator, and it was shown in Refs. |I|J27|J28[| that 
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expression ( P5a|) remains essentially valid if the chaotic part of the emission function is 
restricted to the “core” of direct pions and short-lived resonance decays, we expect Eqs. 
to be practically useful even when resonance decays are included. 


IV. CONCLUSIONS 

We have studied the question to what extent three-pion Bose-Einstein correlations can 
provide independent information about the space-time structure of the emitting source which 
cannot be extracted from two-pion correlations. For chaotic sources we found that the three- 
pion correlator depends on the phase of the two-particle exchange amplitude which drops 
out from the two-particle cross section. This phase can be isolated by proper normalization 
of the true three-pion correlator with respect to the two-pion correlator. It was shown to be 
sensitive to the momentum dependence of the point of highest emissivity in the source and 
to the asymmetries of the emission function around that point. However, this sensitivity is 
weak (it enters only at 4th order in the relative momenta g^), and the corresponding source 
properties are hard to measure. 

We then proceeded to study sources which are not completely chaotic but contain a 
coherent component. We showed that in this case the emission function can be written as a 
sum of two Wigner densities describing the chaotic and coherent components, respectively, 
and expressed the two- and three-pion correlation functions via these chaotic and coherent 
Wigner densities. We showed that a comparison of two- and three-pion correlators allows for 
a determination of the degree of coherence in the source, without contaminations from reso¬ 
nance decays. To this end one must study the respective correlation functions at vanishing 
relative momenta of all particles. To facilitate the extraction of this limit from experimental 
data we derived in Eqs. (|3^ ) and (|3^ ) the most general parametrizations for the two-and 
three-pion correlation functions at small relative momenta. These new parametrizations are 
based on our expressions of the correlation functions in terms of the Wigner density of the 
source; they are exact up to second order in the relative momenta, i.e. for emission functions 
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S{x, K) with a Gaussian x-dependence. After eliminating the redundant g-components, they 
are seen to depend on 16 parameters which are all functions of the average momentum K 
of the pion pair resp. triplet. To determine all these parameter functions, a complete study 
of the two-and three-particle spectra as functions of all 6 -|- 9 = 15 momentum components 
is necessary. (The 16th parameter, e{K), describes the degree of coherence and enters the 
normalization of the correlation functions at vanishing relative momenta.) This is certainly 
not an easy task, and it might be worthwhile to study whether, for certain simple but 
not too unrealistic models for the emission function, it is not possible to obtain simpler 
parametrizations (for example by exploiting certain symmetries of the source). 

Our results show that in the case of partially coherent sources the three-pion correlator 
contains independent information on the second space-time moments of the source which 
cannot be extracted from the two-pion correlator. This information is needed to separate 
the space-time characteristics (lengths of homogeneity or effective widths) of the chaotic and 
coherent parts of the emission function. To extract it in practice will not be easy, but the 
theoretical framework by which this should be done has been presented here. 
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